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ABSTRACT 

This  thesis  analyzes  the  performance  of  a  ramjet  with  an 
Oswatisch  inlet  using  a  blunt  centerbody  and  compares  perfor- 
mance to  a  baseline   ramjet  using  an  inlet  with  a  conical 
spike  at  Mach  3.0. 

Inlet  performance  as  a  ratio  of  inlet  lip  to  nose  center- 
body  ratio,  rT /r  ,  is  developed.   The  capture  streamline  for 
each  ratio  is  determined  and  the  coefficient  of  additive  drag 
is  calculated  as  a  function  of  rT /r  .   Setting  thrust  coef- 
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ficient  equal  to  coefficient  of  drag,  the  performance  of  two 
ramjets  is  determined.   One  ramjet  is  the  baseline  with  a 
spike  inlet;  the  other  ramjet  uses  the  blunt  centerbody. 

Ramjets  and  inlets  are  compared  on  the  basis  of  specific 
fuel  consumption,  excess  thrust  coefficient  and  specific 
thrust.   For  the  ramjet  with  blunt  centerbody,  performance 
parameters  were  calculated  as  a  function  of  inlet  lip  radius 
to  nose  centerbody  radius.   Also  compared  is  the  effect  of 
the  ratio,  r»/r  ,  on  relative  detection  range.   For  both 
types  of  ramjets,  the  detection  range  is  reduced  by  approxi- 
mately 66%.   Performance  of  the  ramjet  with  blunt  nosed 
centerbody  is  severly  handicapped  due  to  high  additive 
drag  and  poor  pressure  recovery.   Specific  fuel  consumption 
is  approximately  50%  greater  for  the  ramjet  with  the  blunt 
centerbody  compared  to  the  ramjet  with  the  spike  inlet. 
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I.   INTRODUCTION 

A.   BACKGROUND 

The  rapid  development  and  operational  deployment  of  long 
range  antiship  cruise  missiles  has  generated  demanding 
requirements  for  the  anti-ship  missile  defense  (ASMD)  mis- 
sion.  The  threat  has  necessitated  that  the  Navy  develop  the 
concept  of  a  "defense  in  depth"  for  the  defensive  doctrine 
of  a  carrier  battle  group  (CBG) .   The  doctrine  emphasizes 
the  utilization  of  all  available  weapon  systems  in  a  layered 
defense  to  defeat  hostile  targets.   Conceivably,  an  anti- 
ship  missile  (ASM)  could  survive  attacks  from  extremely  long 
range  missiles,  the  E-2C  and  F-14  combat  air  patrol  (CAP) 
team,  Standard  ER  and  Standard  MR  missiles,  conventional  gun 
ordnance,  point  defense  weapons  and  close  in  defense  weapons, 
such  as  Phalanx,  prior  to  reaching  a  high  value  unit.   Since 
enemy  doctrine  is  to  fire  the  ASM  in  such  numbers  as  to 
overwhelm  the  Command  and  Control  capability  of  the  CBG 
and  of  an  individual  vessel,  the  accent  of  new  ASMD  weapons 
systems  is  on  quick  response  time,  short  time  of  flight  and 
a  high  degree  of  accuracy  against  a  manuvering  ASM. 

Although  a  single  saturation  raid  could  possibly  be 
defeated,  what  about  a  second  or  third  raid  in  a  single  day? 
Will  the  CBG  have  enough  missiles  to  supplement  airborne  CAP 
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and  still  destroy  the  ASM  threat  30  to  40  miles  from  the 
high  value  unit? 

One  solution  to  this  vexing  problem  is  to  modify  the 
5"/54  Semi-Active  Laser  Guided  Projectile  (SALGP)  which  is 
illustrated  in  Figure  1-1.   The  SALGP  round  is  currently  un- 
dergoing operational  evaluation  (OPEVAL)  prior  to  fleet  wide 
introduction  [Ref.  1]  and  consists  of  a  five  inch  projectile 
with  a  laser  seeker  on  the  projectile  nose.   The  round  is 
designed  to  "home  in"  on  laser  energy  reflected  from  a  target 
by  a  laser  illuminator.   The  modification  to  the  SALGP  round 
consists  of  replacing  the  laser  seeker  head  with  an  infrared 
(IR)  seeker  and  adding  a  solid  fuel  ramjet  to  power  the  pro- 
jectile.  The  ramjet  would  provide  constant  speed  over  a 
longer  portion  of  the  flight  path  and  provide  energy  in  the 
terminal  phase  of  flight  to  allow  for  manuvering  to  hit  a 
manuvering  target.   The  modified  SALGP  is  identified  as  a 
gun  launched  missile  (GLM) .   The  modifications  would  allow 
the  GLM  to  become  a  fire  and  forget  round  requiring  no  ter- 
minal guidance  from  the  ship. 

Although  not  helping  in  the  long  range  defense  of  a  CBG, 
the  GLM  would  complement  the  point  defense  and  close  in 
defense  weapon  systems  to  the  ships  in  the  U.S.  Navy  today. 

The  GLM  will  allow  a  naval  vessel  without  a  missile  sys- 
tem and  with  a  properly  configured  gun  system  to  increase 
its  anti-air  warfare  (AAW)  capability  in  a  cost  effective 
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way.   More  GLM  rounds  can  be  carried  in  existing  gun  maga- 
zines than  missile  rounds  can  be  carried  in  current  missile 
magazines.   Within  a  ten  nautical  mile  zone  around  a  ship, 
the  GLM  will  have  a  faster  response  time,  a  shorter  time  of 
flight  and  be  as  accurate  as  guided  missiles  in  the  current 
Navy  inventory.   Fast  response  and  short  time  of  flight  are 
due  to  the  initial  velocity  imparted  to  the  GLM  by  the  5"/ 54 
gun.   The  shape  of  the  nose  of  the  GLM  round  will  be  con- 
strained by  the  use  of  the  hemispherical  lense  surface  for 
the  infrared  seeker  in  the  projectile  nose.   This  blunt  nose, 
although  optimized  for  the  sensor  optics,  is  not  necessarily 
the  best  aerodynamic  shape. 

Although  the  SALGP  has  the  potential  to  be  highly  effec- 
tive in  the  Naval  Gunfire  Support  (NGFS)  and  in  the  Surface 
to  Surface  Warfare  (SUW)  missions,  it  was  not  designed  ini- 
tially for  the  Air  Warfare  (AAW)  or  Anti-ship  Missile 
Defense  (ASMD)  missions.   Consider  a  GLM  based  on  SALGP 
hardware  to  the  maximum  extent  possible.   To  become  success- 
ful in  the  AAW/ASMD  mission  areas,  it  is  essential  to  im- 
prove the  overall  SALGP  performance  and  capability.   One  area 
to  achieve  easily  measured  performance  gains  is  in  the 
method  of  propelling  the  SALGP. 

One  propulsion  proposal  that  decreases  mission  flight 
time  and  increases  mission  effectiveness  is  to  replace  the 
solid  propellant  rocket  on  the  SALGP  round  with  that  of  a 
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ramjet  with  either  liquid  or  solid  fuel.   Obviously,  to  keep 
total  system  developement  and  procurement  costs  low,  and 
therefore  warrant  the  use  of  a  ramjet  in  a  cost  effective 
analysis,  the  SALGP  round  should  be  changed  as  little  as 
possible  in  adding  the  ramjet.   Effectively,  the  designer  is 
constrained  to  use  a  blunt  nosed  centerbody  to  house  the  IR 
seeker  and  sensor  package.   The  remainder  of  the  projectile 
is  constrained  by  the  dimensions  of  the  handling  equipment 
in  the  MK  45  mount  and  Department  of  the  Navy  Specifications 
[Ref.  2]. 

B .   PURPOSE 

One  purpose  of  this  thesis  is  to  investigate  the  effect 
that  a  blunt  nosed  centerbody  has  on  the  pressure  recovery 
of  an  inlet.   The  pressure  recovery  is  an  important  indica- 
tor of  ramjet  performance,  for  poor  pressure  recovery  will 
reduce  the  flight  envelope  over  which  the  ramjet  operates. 

Although  drag  is  large,  the  blunt  nosed  body  was  aerody- 
namically  feasible  for  a  rocket  guided  projectile,  since  the 
operation  of  the  rocket  engine  does  not  require  air  to  be 
"brought  onboard*1.   The  use  of  a  ramjet  may  prove  inadvise- 
able  or  infeasible  due  to  aerodynamics  of  the  flow  around 
the  blunt  nose  body  enroute  to  the  inlet  annulus .   The  pres- 
sure recovery,  or  lack  of,  may  also  require  the  body  shape 
surrounding  the  IR  sensor  system  be  reconfigured,  within 
existing  system  constraints,  to  allow  for  ramjet  operation. 
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The  aerodynamic  model  used  was  that  of  a  hemispherical  body 

with  an  attached  cylindrical  body  of  constant  diameter.   The 

flow  was  calculated  for  a  cylinder  length  equal  to  four 

times  the  nose  radius.   Distances  involved  were  normalized 

to  the  nose  radius,  r  .   This  is  the  same  as  the  radius  of 

n 

the  hemisphere. 

The  National  Aeronautics  and  Space  Administration  Ames 
Research  Center  (NASA  Ames)  calculated  the  flow  over  a  hem- 
ispherical body  with  an  attached  cylinder  over  the  range  of 
Mach  numbers  from  1.8  to  3.4,  in  steps  of  0.2,  at  zero  angle 
of  attack  and  at  Mach  3.4  with  an  angle  of  attack  (a)  equal 
to  10  degrees.   The  calculations  provide  flow  field  data 
which  included  shock  shape,  sonic  line,  velocity  components, 
pressure,  density,  entropy  and  similar  flow  parameters.   The 
calculations  were  done  by  computer  code  contained  in  two 
separate  programs.   For  a  further  discussion  of  the  programs 
and  the  various  inputs  see  Appendix  A. 

Section  II  reports  on  inlet  performance  and  the  inlet 
flow  field.   Section  III  reports  on  ramjet  performance  for  a 
blunt  nosed  centerbody  for  the  inlet  and  compares  the  per- 
formance of  a  ramjet  with  an  isentropic  spike  inlet. 
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II.   FORMULATION  OF  INLET  PERFORMANCE  EQUATIONS 
AND  INLET  FLOW  FIELD 


A.   INLET  DESIGN  OPTIONS 

In  developing  a  quantitative  basis  for  comparison  of  var- 
ious inlet  designs,  several  items  of  information  are  needed. 

First,  information  needed  is  the  drag  at  the  different 
Mach  numbers  and  a  correlation  of  drag  versus  Mach  number. 
The  blunt  nosed  centerbody  produces  a  detached  shock  that 
stands  off  a  distance  from  the  body.   A  conical  inlet  pro- 
duces an  attached  oblique  shock. 

Second,  the  on-design  performance  of  the  inlet  design 
must  be  evaluated.   In  particular,  the  pressure  recovery  of 
the  inlet  is  important.   What  happens  to  this  recovery  ratio 
when  the  position  of  the  annulus  is  moved  axially  along  the 
body?   If  the  location  is  fixed  and  the  radius  of  the  inlet 
lip  varied  what  happens  to  the  inlet  performance?   A  compar- 
ison must  be  made  of  the  relative  ranges  attained  by  a  ram- 
jet with  an  isentropic  spike  and  one  with  a  blunt  nosed 
centerbody  must  be  made.   Also  of  concern  is  how  sensitive 
the  ramjet  performance  is  to  angle  of  attack  variations. 

Third,  off -design  performance  of  the  ramjet  is  of  inter- 
est in  several  aspects.   First,  how  does  the  pressure  recovery 
of  the  inlet  vary  at  different  Mach  numbers?   Ramjet  perfor- 
mance is  sensitive  to  the  inlet  pressure  recovery.   Generally 
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ramjet  performance  tends  to  improve  initially  at  higher  Mach 
numbers.   Second,  how  does  this  affect  projectile  drag  and 
what  effect,  if  any,  is  transmitted  to  the  projectile  range? 
Third,  what  is  the  operational  region  over  which  the  inlet 
and  ramjet  will  continue  to  operate?   The  lower  performance 
limit  of  the  inlet  is  that  Mach  number  where  the  normal  shock 
within  the  inlet  is  expelled.   This  condition  is  refered  to 
as  subcritical  ramjet  operation  and  can  lead  to  a  condition 
called  "buzz".   The  upper  performance  limit  of  the  inlet  is 
the  maximum  flight  speed  at  which  adequate  pressure  recovery 
can  be  attained.   Last,  how  does  the  angle  of  attack  affect 
off  design  performance? 

The  desire  that  the  ramjet  powered  GLM  fit  within  the 
existing  confines  of  the  5"/54  MK  45  gun  limits  the  inlet 
design  to  that  of  the  axisymmetric  variety.   In  theory, 
"pop-up"  inlets  or  scoops  could  be  used.   However,  the  ability 
of  the  inlets  to  survive  a  gun  launch  of  8000  g's  and  func- 
tion reliably  increases  the  complexity  of  the  round;  variable 
geometry  inlets  are  not  a  good  use  of  the  remaining  available 
space  in  the  projectile.   The  small  diameter  of  the  projectile, 
coupled  with  the  nose  mounted  laser  seeker,  requires  intel- 
ligent use  of  available  space. 

The  axisymmetric  inlets  which  could  be  utilized  are 
1)  blunt  nosed  centerbody,  2)  conical  nose  and  variations, 
3)  an  unconstrained  design,  an  isentropic  spike  for  example, 
and  4)  a  cylindrical  analog  of  the  Busemann  Biplane. 
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The  blunt  nose  shown  in  Figure  II -1  is  currently  in  use 
on  the  SALGP  projectile,  without  the  annul us .   As  modeled  by 
the  simulation,  the  inlet  consists  of  a  hemispherical  nose 
connected  to  a  constant  diameter  body  extending  beyond  the 
inlet  lip.   The  conical  nose  in  Figure  1 1  -  2  uses  a  cone  of 
fixed  angle  to  focus  the  shock  wave  at  the  inlet  lip  when 
at  the  design  Mach  number.   A  variation  is  shown  in  Figure 
1 1 -  3 .   Shown  in  Figure  1 1 -4  is  the  unconstrained  design. 
The  design  would  consist  of  an  isentropic  spike.   Figure 
1 1 - 5  shows  the  cylindrical  analog  of  the  Busemann  Biplane. 
The  inlet  "swallows"  the  oblique  shock  wave  formed  by  the 
nose  at  the  design  Mach  number  and  has  a  body  of  constant 
diameter.   The  advantage  is  the  elimination  of  external  wave 
drag . 

B.   CALCULATION  OF  STAGNATION  PRESSURE  RATIO  OF  THE  INLET 

For  calculation  purposes,  the  flow  around  the  blunt  nosed 
body  is  divided  into  two  segments  as  shown  in  Figure  I I -6. 
Region  1  is  the  region  of  flow  from  upstream  infinity  to  the 
shoulder  of  the  centerbody.   The  centerbody  shoulder  is 
located  in  plane  S  in  Figure  1 1  -  6 .   Region  1  is  characterized 
by  flow  in  both  the  longitudinal  and  radial  directions. 
Region  2   is  the  region  of  flow  downstream  of  plane  S.   The 
majority  of  the  flow  is  in  the  longitudinal  direction  with  a 
relatively  small  radial  component.   The  origin  for  coordi- 
nate systems  for  regions  1  and  2   is  the  same  as  illustrated 
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Figure  II-l   Supersonic  Inlet  With  Blunt  Centerbody 


Figure  1 1 -  2   Supersonic  Inlet  With  Conical  Spike 
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Figure  1 1 - 3   Supersonic  Inlet  (variation) 


Figure  1 1 - 4   Supersonic  Inlet  With  Isentropic  Spike 
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Figure  1 1 - 5   Cylindrical  Analog  Of  Busemann  Biplane  Inlet 
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Figure  I I- 6   Flow  Regions  Around  A  Blunt  Nosed  Body 
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in  Figure  1 1  -  6 .   The  coordinates  for  region  1  are  X  and  Y, 
which  are  the  symbols  used  in  the  NASA  Ames  computer  printout. 
The  coordinates  for  region  2  are  R  and  Z,  also  illustrated 
in  Figure  II-6.   Each  of  the  regions  has  different  reference 
values  and  symbols  for  the  flow  properties.   Appendix  A  lists 
the  reference  values  and  symbols. 

In  a  real  ramjet,  the  sizing  of  the  inlet  capture  area, 
throat  area  and  inlet  area  can  dramatically  vary  the  perfor- 
mance at  different  Mach  numbers.   The  importance  of  the  effect 
the  inlet  plays  on  performance  is  paramount. 

Figure  1 1 - 7  shows  the  magnitude  of  pressure  losses  as  a 
fraction  of  total  pressure  recovery.   As  can  be  clearly  seen, 
the  majority  of  loss  is  due  to  boundary  layer  growth.   How- 
ever, as  the  flight  Mach  number  increases,  note  that  all  the 
losses  increase  in  magnitude. 

For  a  conical  inlet,  the  oblique  shock  losses  can  readily 
be  calculated.   The  overall  pressure  recovery  value  of  0.72, 
for  Mach  3.0  flow,  can  be  obtained  from  Figure  II-7. 

For  the  blunt  nosed  body,  p./p    is  tabulated  in  the  NASA 
Ames  printout.   This  corresponds  to  Pt^/Pri'  t^ie  °bllclue 
shock  losses,  in  Figure  1 1  -  8 .   To  account  for  the  varying 
pressure  recovery  values  within  the  annulus ,  a  mass  weighted 
average  will  be  used.   This  concept  will  be  discussed  in 
Section  II-C-3. 
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Magnitude  Of  Pressure  Losses 
(Reproduced  from  Ref.  4,  p.  64) 
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Total  pressure  ratio,  or  pressure  recovery,  is  defined 
as  77 , .   From  Figure  1 1  -  8  , 


PtS     H2     ^tS  n, 

"d  =  P^  =  Ptl  Pt2  tl} 


By  defining  tt,'  as  equal  to  Pt.?/Pti»  tne  Pressure  recovery 
becomes 


pt5 
*d   =   V  p~  <2> 


The  value  of  Pte/P*2  can  ^e  determined  from  Figure  1 1  -  7 . 
The  value  of  Ptc/Pt2  represents  the  accumulation  of  losses 
in  an  inlet  due  to  the  subsonic  diffuser  loss,  Ptc/PtA> 
transonic  shock  and  viscous  losses,  P^/p^tJ  and  pressure 
loss  due  to  boundary  layer  growth,  p  .,/p  - .   The  value  of 
Ptq/Pt-2  *s  °b'taine<i  from  Figure  II -7  as  a  ratio  of  overall 
inlet  pressure  recovery  to  oblique  shock  losses.   For  Mach 
3.0  flow 

PtS     0  72 

-i2  =  ^4t  =  0.7423 

pt2      0.97 

The  value  of  ir,'  was  calculated  using  the  flow  field  calcu- 
lations in  the  NASA  Ames  printout.   The  overall  effect  of  the 
reduced  recovery  ratio  is  that  the  engine  efficiency  and  the 
pressure  ratio  at  the  nozzle  exit  is  changed.   In  turn,  the 


31 


change  in  the  pressure  ratio  may  cause  over  or  under  expanded 
flow  from  the  exit  nozzle  [Ref .  3] . 

C.   FLOW  FIELD  DESCRIPTION  IN  THE  VICINITY  OF  A  BLUNT  NOSE 
For  the  blunt  nosed  centerbody  at  zero  angle  of  attack, 
the  flow  field  is  axisymmetric  with  respect  to  the  centerline 
of  the  projectile.   The  calculation  of  pressure  recovery  and 
ramjet  performance  is  greatly  simplified. 
1 .   Determination  Of  Capture  Streamline 

Inherent  to  the  analysis  of  ramjet  performance  is  the 
determination  of  the  mass  flow  rate  of  air  at  the  annulus  of 
the  inlet,  nu  .   After  determining  the  mass  flow  rate,  m. , 
the  streamline  coordinates  will  be  determined  at  various 
points  by  equating  the  local  mass  flow  to  the  mass  flow  at 
the  lip.   From  these  coordinates,  the  additive  drag  and 
coefficient  of  additive  drag  for  the  projectile  at  a  parti- 
cular Mach  number  will  be  determined.   Two  approaches  for 
calculating  the  streamtube  points  were  tried  and  will  be 
discussed  here. 

a.   Mass  Flow  Method 

The  first  method  involves  mass  flow  analysis. 
The  mass  flow  through  any  cross  section  of  the  particular 
streamtube  will  equal  that  of  the  mass  flow  through  the 
annulus.   Of  particular  interest  is  the  radius  of  the  capture 
area,  r  ,  on  the  bow  shock  wave.   The  air  captured  by  this 
area  goes  into  the  inlet.   The  remainder  of  the  air  spills 
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past  the  inlet  and  affects  the  performance  of  the  projectile 
by  adding  to  the  total  vehicle  drag  with  a  term  called  addi- 
tive drag,  D  . 

a 

The  mass  flow  at  the  lip  is  given  as 


mL   =  J   PL   u 


mT   =   /  pT   u.  dA  (3) 


Where  py  is  the  density,  dA  the  elemental  area  and  uT  the 
velocity  component  normal  to  that  area.   From  the  continuity 
equation,  the  mass  flow  at  the  shock  front  is 


2 
m   =    p  iTr  u  (4) 

s       s   s   s 


Refer  to  Figure  1 1 - 9 .   By  equating  equations  3  and  4 

lip 

PL 
body 


/ 


P   Trr      u  =        /    pT      uT    27TrT    drT  (5) 

s      s      s  ILL  L        L 


By  normalizing  the  distances  involved  to  the  nose  radius  and 
by  dividing  by  p^  and  v  ,  to  allow  use  of  the  NASA  Ames 
printout  values  in  region  2  ,  the  equations  become 


"2 =   2  p—    /p2  u2  r   dr  (6) 

r  v  M«>    J.     , 

n  m  body 
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where  u2   =  uL/vm;   P2   =   PL/Ptoo;  and  r  =  r/i^.   The  primed 
values  are  the  tabulated  values  from  the  NASA  Ames  printout 
and  the  subscript  refers  to  the  applicable  flow  region  de- 
fined in  Figure  II -6.   Q   is  defined  as  u  /v  ,  and  equation  6 

a  ^oo  S    m 

becomes 


2       -  ,  . 

rs          2  pt~  r1^    ' 

T  =   Qj-F  /  P2   U2 

n  body 


r   dr  (7) 


The  integrand  in  equation  7  was  plotted  and  the  area  under 
the  curve  calculated  for  the  desired  lip  radius.   Since 
r   =  1.0  inches,  Q  =  0.80174,  p  /p^   =  0.07623  at  Mach  3  and 
the  integral  from  body  to  lip  is  0.00805,  r   can  be  determined 
For  Mach  3  airflow,  with  rT  =1.2 


2  2 

rs   =    (0.80174) (0. 07623)   CO. 00805) 


Hence,  r   =  0.5133. 
s 

Table  I I - 1  shows  r   for  Mach  3.0  and  a  nose  centerbody  radii 

s  ' 

equal  to  1.00  inches.   The  maximum  value  of  the  lip  radius 
is  1.854  inches.   This  is  the  size  imposed  by  the  constraints 
of  the  5"/54  handling  system.   Now  the  streamline  locations 
at  different  flow  field  points  can  be  calculated.   The  compu- 
ter program  used  is  described  in  Appendix  B.   The  flow  geometry 
is  shown  in  Figure  I I -9. 
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Table  II-l 

Capture  Radius,  r  ,  As  A  Function  Of  Lip  Radius,  i\  ,  For  r 

Equal  To  1.0  Inch  At  Station  Z=1.00 

rT  (inches)  r   (inches) 

1.1000  0.3821 

1.2000  0.5883 

1.3000  0.7765 

1.4000  0.9610 

1.5000  1.1474 

1.6000  1.3353 

1.7000  1.5231 

1.8000  1.7066 

1.8540  1.8086 
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The  mass  flow  rate  across  a  line  of  constant 
angle,  a,  is  given  by  the  following  equation 


/cs 


ma  =   /  p  u„  dS  (8) 


where  b  refers  to  the  body  and  cs  refers  to  the  capture  stream- 
line in  Figure  1 1  -  9 .   From  Figure  11-10,  an  element  of  area 
on  the  line  of  constant  a  is  given  as 

dS   =   27Trdl  (9) 

and  dl  can  be  related  to  the  angle  a  by 


dl    =   -£l_  (10) 

sma  *■  J 


Hence,    dS    in   terms   of  dr   becomes 


dS        =        27rr  -^—  (11) 

sma  *■      J 


From  Figure  11-11,  an  angle  3  is  defined  so  that 


tan_1(v/u)  (12) 
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X 


Figure  11-10  Magnified  View  Of  Geometry  Between 
Points  1  And  2  Of  Figure  II -9 
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Figure  11-11   Geometry  To  Determine  Component  Of  Velocity 
Normal  To  The  Area  For  Integration  Of  Mass 
Flux 
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Hence,  the  normal  component  of  the  velocity  becomes 


/    2    2 
u    =  /   u  +v   cos  (ir/2-a-  3)  (13) 

n 

With  the  substitution  of  a  trigonometric  identity  for  the 
cosine  term,  equation  13  can  be  simplified  to 

I 2 2~~ 

u    =  /   u  +  v   sin(ct+  3)  (14) 

The  geometry  of  Figure  11-11  gives 

cos3   =       u  (15) 

/  u  +v 

and 


sin  3   =       v  (16) 

/  u  +v 


Using  equations  15  and  16  and  substituting  into  equation  14 
leaves 


u    =   /  u2+v2   (u  sina  *   v  cosa   )  (17) 

/  u  +v     /  u  +v 


0 
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The  value  of  airflow  normal  to  the  cross  sectional  area  dS 
then  becomes 


u    =   u  sina  +  v  cosa  (18) 

n 


Substituting  equations  11  and  18  into  equation  8,  the  follow- 
ing equation  is  obtained 


m    =   /  (p2irr)  (u  sina  +  v  cosa)  /(sina)  dr 


This  becomes ,  after  manipulation 

cs 
m  =     /     p    2-rrr      (u      +   v      cota)dr  (20) 


m    =  I    p  2tjt   (u  +  v  cota^ 


The  mass  flow  given  by  the  above  equation  through  the 
streamtube  at  the  radial  line  of  constant  angle  a  must  equal 

the  mass  flow  through  the  nose  capture  area.   Multiplying 

2 
the  integrand  by  (r  /r  )   and  dividing  both  sides  by  the  mass 

flow  at  the  bow  shock,  equation  4,  gives  the  following  result 


■/; 


cs      r 

2    (-J!)2   £-     L_      [U_  +   v_]    ll                               (21) 

vrprLu  ur 

s          s        n          s  s        n 
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In  front  of  the  projectile  shoulder,  equation  21  becomes 

cs 

(22) 


(-£-)     =    2   /  p±   r   (ux  +  v1  cota)dr 
n  J  x. 


For  region  2  ,  which  is  after  the  shoulder,  due  to  the  dif- 
ferent reference  values  equation  20  becomes 

2p 


s,2 


t°°     c    •    *  ^  •       • 

—   /   p~  r   (u-  +  v-  cotajdr 
P„   7    2     ^2    2 


Cr2")     =   ?|-S-   /   P?  r   (u^  +  v9  cota)dr       (23) 

"b 


where  in  equations  22  and  23,  the  primed  values  are  taken 
from  the  NASA  Ames  printout.   In  region  2  ,  a  equals  90°  and 
cotangent  of  a  is  zero. 

b.   Angular  Method 

The  second  method  of  streamtube  calculation  is 
by  the  angular  method.   A  linear  variation  of  the  known  flow 
field  properties  over  small  distances  is  assumed.   The  method 
starts  at  the  lip  and  progressively  moves  forward  along  the 
body  determining  streamtube  position. 

Initially,  the  local  velocity  vector  at  the 
various  flow  field  points  must  be  plotted  in  polar  notation. 
By  averaging  the  velocity  flow  angles  at  two  adjacent  stations, 
an  average  angular  value  can  be  determined.   As  an  example, 
for  a  lip  radius  of  1.2 
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934    =    (93  +  94)/2  (24) 


where   0,  =  local  flow  field  velocity  angle  at  station  3 

0.  =  local  flow  field  velocity  angle  at  station  4 

0_4=  average  angle  between  station  3  and  4  stream 
tube  locations 

9.,  the  flow  field  angle  at  the  lip,  is  determined 

from  linear  interpolation  of  the  known  flow  field  angles.   By 

varying  9.,,  with  9.  known,  the  value  of  9,  is  determined. 

As  a  function  of  9.,. 


94   =    2943  '  93  (25) 


9.,  is  then  plotted  as  a  function  of  the  radius  of  the 
streamtube,  R^  in  this  case. 

By  studying  Figure  11-12,  the  geometry  of  the 
field  also  determined  9.,. 

R.  -  R- 
tan  943   =   [  z4  .  z   ]  (26) 

In  the  above  equation,  all  quantities  with  the  exception  of 
R3  are  known.   9.,  is  then  plotted  versus  R.,,. 
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Figure  11-12   Flow  Field  Geometry  For  The  Angular 
Method  Of  Streamtube  Determination 
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The  intersection  of  these  two  curves  determines 
R^  and  9*7,  and  hence  9^  can  be  determined  by  interpolation. 
In  this  manner,  the  streamtube  points  are  calculated  for  the 
desired  lip  radius.   This  method  produces  comparable  results 
to  the  mass  flow  analysis  but  is  not  accurate  enough  pro- 
ceeding past  the  nose  of  the  projectile.   Additionally,  the 
process  is  extremely  slow  and  complicated. 
2 .   Determination  Of  Additive  Drag 

The  additive  drag  is  calculated  by  a  computer  program 
discussed  in  Appendix  C.   Assuming  that  the  points  defining 
the  edge  of  the  streamtube  determined  by  mass  flow  analysis 
can  be  connected  by  straight  line  segments,  the  equation  for 
additive  drag  can  be  written  as 


D„    =    Z(p.  -  p  )A.  sin0.  (27) 

a        Vii   ^°°J    li  K      J 


where  p.  equals  the  average  pressure  at  the  data  points  i 
and  i+1;  0.  is  the  angle  of  the  straight  line  segment  between 
the  points  i  and  i+1  measured  from  the  projectile  axis;  and 
A.  is  the  surface  area  of  the  cone  between  the  i  and  i+1  data 
points.   The  geometry  is  shown  in  Figure  11-13. 

The  coefficient  of  additive  drag  (C,  ,)  can  be  de- 
termined by  dividing  equation  27  by  q^  and  reference  area,  A  , 

pi  r 

2  2  pM(^i  -  1)  A.  sin0. 

Cdad   =   -  2 (28) 

YP   M   A 

1  r  00     oo       J* 
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Figure  11-13   Geometry  For  Additive  Drag  Calculation 
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The  reference  area,  A  ,  is  taken  to  be  the  base  of  the  five 
inch  projectile.   Simplification  of  equation  (2  8)  leads  to 


Cdad   "   -TT7  E  (F  "  1}  Sr  sin9i  W 

yM        °°       r 


Equation  (29)  is  the  equation  programmed  in  Appendix  C. 

The  program  utilizes  the  points  generated  by  the  pro- 
gram in  Appendix  B  to  calculate  C,  ,.   However,  it  is  impor- 
tant to  scale  the  nose  capture  radius  from  a  centerbody  radius 
of  one  inch  to  the  actual  centerbody  radius  on  the  projectile. 
Previously  rT  was  fixed  at  1.854  inches.   A  ratio  between 

the  centerbody  and  lip  radii  is  defined  as  x  =  rT/r  . 

Li      n 


The 


lip  radius,  r, ,  is  fixed  by  the  dimensions  of  the 


handling  equipment,  and  r  can  be  determined  as  a  function  of 

the  ratio  x.   Table  1 1 - 2  displays  the  values  of  x,  r  ,  C,  j, 

and  mass  flow  ratio.   Mass  flow  ratio  is  defined  as  m  /m  . 

s   s 

The  symbols  m  and  m  are  defined  in  Table  II-2. 

3 .   Average  Stagnation  Pressure  Ratio  At  The  Inlet  Lip 
The  stagnation  pressure  at  the  annulus  of  the  inlet 
is  a  function  of  radius.   A  mass  weighted  average  stagnation 
pressure  was  calculated.   The  average  pressure  recovery  will 
be  multiplied  by  a  factor  of  ptr/pt2  t0  obtain  total  pressure 
recovery  as  discussed  in  section  II. 
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Table  II-2 

Summary  Of  Radius  Of  Capture  Streamtube,  Additive  Drag  Coef 

ficient,  And  Mass  Flow  Ratio 


For  r.  ■  1.854  in. 


x=rT  /r 
L  n 


r  /r 
s   n 


n 


'dad 


m  /m 
s   s 


1.1 

0.3821 

1.6855 

1.2 

0.5883 

1.5450 

1.3 

0.7765 

1.4262 

1.4 

0.9610 

1.3243 

1.5 

1.1474 

1.2360 

1.6 

1.3353 

1.1588 

1.7 

1.5231 

1.0906 

1.8 

1.7066 

1.0300 

1.854 

1.8086 

1.0000 

0.6440 
0.9089 
1.1074 
1.2727 
1.4181 
1.5473 
1.6611 
1.7578 
1.8086 


0.3726 
0.2956 
0.2299 
0.1775 
0.1319 
0.0928 
0.0588 
0.0304 
0.0179 


0.121 
0.240 
0.357 
0.471 
0.595 
0.697 
0.803 
0.899 
0.952 


m  is  defined  by  equation  (4) ;  m°  is  obtained  from  equation 


(4)  with  r  replaced  by  rT 
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The  mass  weighted  average  of  pressure  recovery  is 
given  by 

Pt  x        Pt  mi 

tt,      =      —      =      s     (— ) .  r1  (30) 

Ft°°  j=l      rt°°   J      L 

where  the  subscript  j  refers  to  the  jth  streamtube  within 
the  inlet,  m.  is  the  mass  flow  of  the  jth  streamtube  and  mT 
is  the  total  mass  flow  through  the  entire  inlet  as  defined 
by  equation  3.   This  geometry  is  further  shown  in  Figure 
11-14.   If  the  inlet  lip  splits  a  streamtube  the  values  for 
pressure  recovery  and  mass  flow  are  linearily  interpolated 
to  achieve  the  average  pressure  recovery.   The  values  of 
(p  /p   ) .  were  taken  from  the  NASA  Ames  printout. 
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Outer  Streamline  of  j th  Streamtube 
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Figure  11-14   Geometry  Of  Streamtubes  At  The  Inlet 
Annulus  Used  To  Calculate  Average 
Stagnation  Pressure,  tt  , ' 
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III.   CALCULATION  OF  RAMJET  PERFORMANCE 

A.   PERFORMANCE  EQUATIONS 

Ramjet  engine  performance  can  be  quantified  by  the  specif 
ic   thrust,  specific  fuel  consumption,  and  the  coefficient 
of  thrust.   Performance  equations  will  be  developed  for  each 
of  the  above  performance  criterion  which  include  the  concept 
of  additive  drag. 

If  one  dimensional  flow  is  assumed,  the  thrust,  F,  of  a 
ramjet  engine  is  described  by  Netzer  [Ref.  5]  as 


F  =  mv   -mv   -D  +  A  fp   -  p  )  (31) 

e  e    oo    a    e^e   *o 


where  the  subscript  e  refers  to  exit  conditions  and  the  sub 
script  o  refers  to  the  inlet  conditions.   Term  1  is  defined 
as  jet  or  gross  thrust,  term  2  is  the  ram  drag,  term  3  is 
the  additive  drag  and  term  4  is  the  pressure  thrust.   If 
ideal  expansion  is  assumed  (p   =  p  ) ,  equation  1  becomes 


F=mv   -mv   -D  (32) 

e  e    o  o    a  **  J 


Now,  defining  the  fuel  air  ratio,  f,  as 


f  =  mr/m  (33) 

f  o 
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the  mass  flow  through  the  ramjet  exit  can  be  written  as 


m  =  m  +  m,  =  m   (1  +  f)  (34) 

e    o    f    o 


The  fuel  air  ratio,  f,  is  generally  much  less  than  one. 
Combining  equations  32  and  34  allows  the  total  thrust  to  be 
written  as 


F  =  m  (v   -  v  )  -  D  (35) 

o v  e    oJ  a  *■  J 


Factoring  v  from  the  first  term  on  the  right  hand  side  gives 

v 

F  =  m  v   (-£■  -  1)  -  D  (36) 

o  o  vv     J  a  K      } 

o 


Combining  equation  36  with  the  definition  of  Mach  number  and 
the  speed  of  sound,  a  ,  inserting  an  expression  for  combustor 
energy  balance  and  assuming  that  the  stagnation  pressure 
through  the  inlet,  burner  and  nozzle  are  taken  as  to  be 
constant  allows  thrust  to  be  written  as  a  function  of  fuel 
air  ratio,  heating  value  of  the  fuel  and  the  combustion 
temperature. 


F  =  Vo  [  /  l  +  rr-    -1  1  -  Da  (37) 

p  to 
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The  specific  thrust  is  obtained  from  equation  (37)  by 
dividing  by  the  mass  flow  rate,  m  .   Units  of  m   are  slug/sec 
and  specific  thrust  is  usually  in  units  of  lb,./ (lb  /sec). 


SF   =   t^-  (38) 

mQg 


The  specific  fuel  consumption,  SFC,  is  defined  as  the 
unit  mass  flow  of  fuel  per  hour  per  pound  thrust.   Expressed 
in  the  form  of  an  equation. 


3600mfg 
SFC    =   ^-  (39) 


In  equation  (39)  m.  has  units  of  slug/sec.   Inserting  equation 
37  yields 


3600m.  g 
SFC   =   x (40) 

m  v   [/l  +  3tt^~-  1]  -  D 

o  o  L      C  T_    J     a 
p  to 


where  m-  =  mass  fuel  flow  rate  (slugs/sec) 

g  =  gravitational  constant 

h  =  heating  value  of  fuel  (BTU/lb) 

f  =  fuel  air  ratio 

C  =  specific  heat  capacity  (BTU/lb  °R) 

T  =  inlet  stagnation  temperature 
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The  thrust  coefficient  is  defined  as 


Cf    =    O"  (41} 

^o  r 


where  q   is  the  dynamic  pressure  and  A   is  the  reference 
area.   The  reference  area  is  the  base  area  of  the  five  inch 
projectile.   For  a  ramjet,  mass  flow  at  the  inlet  is  given 
as 


m    =   A  v  o  (42) 

O         0  oo 


Combining  equations  (37)  ,  (41)  and  (42)  gives 

2A„      TT? Do 

C£   =    --   [  /  1  +  t^-    :  -  TT-r- 

The  preceding  equations  were  derived  for  a  ramjet  engine 
without  internal  losses.   Losses  in  a  ramjet  engine  that 
should  be  accounted  for  are 

a)  loss  of  stagnation  pressure  in  the  diffuser,  tt 

b)  loss  of  stagnation  pressure  in  the  burner,  tt, 

c)  loss  of  stagnation  pressure  in  the  nozzle,  tt 

d)  combustion  efficiency,  r\. 


2Ao 

A 
r 

p    to 

•  1] 

D 
a 

q  A 
no   r 

54 


These  terms  result  in  the  generation  of  a  common  term, 


S  ,  which  does  affect  the  calculation  of  ramjet  performance. 


S.  is  defined  as 
r 


Vd^n 


y-1 

)T~  -i 


-*-*—   M 

2    o 


1  + 


n,  fh 
b 

p  to 


1/2 


(44) 


The  square  root  appearing  inequations  (38),  &0)  and  (43)  is  replaced 
by  S   to  give  ramjet  performance  with  internal  losses  and  with 
an  additive  drag  term.   The  performance  equations  becomes 


SF 


[m  M  a   (S   -  1)  -  D  ]/m  g 
L  o  o  o  v  r    '  a   o6 


(45) 


SFC   = 


3600mfg 

— r — 


(46) 


2A  D 

7^   (Sr  "  « 


q  A 
no  r 


(47) 


The  computer  program  developed  by  Fuhs  [Ref.  6]  was 
modified  to  compute  the  additive  drag  term.   Holding  the  values 
of  tt,  ,  TT   and  n,  constant,  the  performance  of  a  ramjet  engine 
with  a  blunt  centerbody  was  predicted.   The  performance  of 
a  ramjet  with  a  conical  spike   inlet  as  predicted  by  the 
program  serves  as  a  baseline  measurement  of  performance. 
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IV.   DISCUSSION  OF  RESULTS 

A.   PERFORMANCE  CRITERION 

Performance  of  a  ramjet  engine  can  be  quantified  in 
different  ways.   For  purposes  of  this  report,  the  performance 

was  compared  to  a  ramjet  baseline  configuration  using  a  coni- 

2 
cal  inlet  spike;  the  inlet  capture  area  is  0.0123  ft   [Ref.  7] 

The  performance  of  a  ramjet  with  the  blunt  centerbody 
configuration  was  calculated  for  a  range  of  fuel/air  ratios. 
Performance  is  specified  on  the  basis  of  combusiton  exit 
temperature  (T.?) ,  specific  fuel  consumption  and  excess 
thrust.   Of  special  interest  is  performance  when  Cr  equals 
C^.   The  value  of  C~,  which  was  obtained  from  White  [Ref.  8] , 
was  0.349. 

The  combustor  exit  temperature  is  calculated  by  the  ram- 
jet performance  program  as  an  output  variable.   Current  tech- 
nology limits  the  steady  combustor  exit  temperature  to  about 
4400°R.   The  value  chosen  for  analysis  purposes  is  4422°R. 
Above  this  temperature,  the  combustor  and  exit  nozzle  will 
melt,  if  run  for  a  continuous  period  of  time.   The  value  of 
T  ,  is  generally  less  than  T  -  at  the  stoichiometric  fuel/air 
ratio  for  the  hydrocarbon. 

For  the  final  seconds  of  flight,  it  is  possible  to  boost 
the  fuel/air  ratio  to  the  stoichiometric  value.   The  boost 
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in  fuel/air  ratio  will  increase  the  amount  of  excess  thrust 
available  for  terminal  maneuvering.   Since  the  GLM  will 
explode  at  flight  termination,  the  combustor  exit  temperature 
is  of  little  importance  in  terminal  flight. 

A  further  basis  for  performance  comparison  is  the  effect 
on  relative  detection  range  of  seeker  aperture.   Seeker 
aperture  is  equal  to  irr 

Specific  fuel  consumption  (SFC)  was  chosen  as  a  compari- 
son standard  because  it  provides  a  good  measure  of  ramjet  fuel 
efficiency.   Excess  thrust  is  used  as  a  measure  of  the  ability 
of  the  guided  projectile  to  maneuver  in  the  terminal  phase  of 
the  encounter.   Regardless  of  the  seeker  accuracy,  lack  of 
thrust  above  the  amount  of  flight  drag  will  result  in  the 
inability  to  pursue  a  maneuvering  target.   Excess  thrust 
coefficient,  Cf  ,  is  calculated  as 


Cfe  (T  3)    =    C£(T  3)  -  CD  (48) 


where  C~  is  airframe  drag  and  C.-CT^)  is  the  thrust  coeffi- 
cient for  T  _.   Two  values  of  T  -  are  of  interest.   One  is 
the  maximum  allowed  combustor  exit  temperature  for  steady 
operation  which  is  4422°R.   The  other  is  T  ,  for  stoichiomet 
ric  combustion  which  is  5968°R;  stoichiometric  fuel/air 
ratio  is  0.06. 
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B.   RESULTS  OF  CALCULATIONS 

Figure  IV-1  plots  T  ,  versus  rT /r  at  a  value  of  C~  =  C,-. 

t  j         l  n  u  i 

As  can  be  clearly  seen,  as  the  ratio  of  rT/r   increases, 

Li  n 

or  as  the  centerbody  shrinks,  T  ,  decreases.   T  ,  reaches 

the  maximum  sustained  operating  value  of  4422°R  when  t*/t 

is   1.63;    rT/r      equal    to    1.63    is   the  minimum  value    for   sus - 
Li     n 

tained  operation.   For  stoichiometric  combustion  (f=0.06), 

rT/r   could  possibly  be  reduced  to  a  value  as  small  as  1.50. 
L  n       c 

From  the  minimum  value  of  rT/r  obtained  in  Figure  IV-1, 

L  n  ° 

Figure  IV-2  is  entered.   Figure  IV-2  is  a  plot  of  excess 

thrust  coefficient  versus  rT/r   and  of  stoichiometric  excess 

thrust  coefficient  versus  rT/r  .   Stoichiometric  excess 

L  n 

thrust  is  defined  as  C,.  at  stoichiometric  conditions  (f=0.060) 
minus  C^. 

Using  a  value  of  rT /r   equal  to  1.63,  the  ramjet  will 

L  n 

have  excess  maneuvering  thrust  only  at  values  of  rr /r   greater 

&  J  L  n  & 

than  1.63.   Use  of  stoichiometric  fuel/air  mixture,  during 

the  last  seconds  of  the  encounter  with  a  maneuvering  target, 

results  in  excess  thrust  being  generated  for  ry /r   greater 

than  1.51. 

Figure  IV-3  is  a  graph  of  the  specific  fuel  comsumption 

(SFC)  versus  rT /r  .   SFC  decreases  as  rT/r   increases.   The 
l  n  l  n 

baseline  ramjet  has  a  SFC  equal  to  2.02  for  C~  equal  to  Cf. 
For  rL/r  equal  to  1.63,  the  ramjet  with  a  blunt  nose  has 
651  higher  SFC. 
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Figure  IV-4  is  a  plot  of  relative  detection  range  of  an 
IR  seeker  for  various  rT /r  .   The  signal-to-noise  ratio  of 
the  IR  seeker  and,  therefore,  the  detection  range  is  propor- 
tional to  the  area  of  the  lense.   The  lense  area  is  directly 

related  to  the  square  of  the  lense  radius.   As  rT/r   in- 

L  n 

2 
creases,  the  detection  range  decreases  as  (r  /r,  )  . 

Mass  flow  ratio  was  defined  in  the  text  near  Table  II-3. 

The  ratio  is  interpreted  as  the  fraction  of  the  mass  flow 

actually  captured,  m  ,  relative  to  mass  flow  into  an  area 

2 
iTr.  .   Figure  IV-5  has  mass  flow  ratio  plotted  as  a  function 

of  rT/r  .   Also  Figure  IV-5  has  a  plot  of  capture  streamtube 

radius  ratio,  r  /rT  ,  as  a  function  of  rT/r  .   For  the  mini- 
'   s   L'  L  n 

mum  value  of  r./r  equal  to  1.63,  the  capture  radius  ratio 
is  0.84.   Likewise,  the  mass  flow  ratio  is  0.70. 

Figure  IV-6  shows  the  variation  of  additive  drag  coeffi- 
cient with  rT/r  .   At  the  minimum  value  of  rT/r   of  1.63, 
L  n  L  n 

C,  ,  is  0.085.   The  airframe  drag  coefficient  is  0.349.   Hence, 

C,  j  is  24%  of  C^.   To  understand  the  behavior  of  C,  A   as  a 
dad  u  dad 

function  of  rT /r   refer  to  Figure  IV- 7 .   For  rT /r   equal  to 
L  n  &  1/  n   n 

1.1,  which  is  illustrated  in  Figure  IV- 7 (a),  the  capture 

streamline  between  bow  shock  and  lip  of  the  annulus  is  steep. 

Also,  the  pressure  is  large  near  the  axis.   Consequently, 

C.j  is  large  for  small  rT /r  . 
dad       b  L  n 

The  greater  the  ratio  of  rT/r  becomes  the  better  ramjet 
&  L  n  J 

performance  becomes.   However,  to  achieve  even  marginal 

performance  at  rT /r   =1.63,  an  unrealistic  penalty  in  the 

l  n 
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a)   rL/rn  =1.1 
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Lip  Of  Annulus 


b)   rT/r   =  1.4 
L'  n 


Figure  IV- 7   Bow  Shock  Wave  And  Capture  Streamtube 

For  rT/r  Equal  To  1 . 1  And  1.4 
L   n   n 
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form  of  a  63%  loss  of  relative  detection  range  must  be  paid. 
Lower  values  of  rT /r   are  unable  to  be  obtained  in  a  thrust 
equals  drag  configuration  due  to  high  values  of  additive 
drag  and  poor  it,. 

Table  IV-1  contrasts  ramjets  with  the  spike  inlet  and 
two  blunt  nose  inlets.   All  three  ramjets  have  identical 
thrust  coefficients  equal  to  airframe  drag  coefficient.   Due 
to  large  specific  thrust,  the  ramjet  with  spike  inlet  re- 
quires considerably  less  mass  flow  rate.   The  ramjet  with 
spike  inlet  has  superior  specific  fuel  consumption.   Excess 
thrust  is  comparable  for  the  three  ramjets. 
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Table  VI-1 
Comparison  Of  Spike  Inlet  Versus  Two  Blunt  Nose  Inlets 


Ramjet  Using 

Parameter 

Spike 
Inlet 

Blunt  Nose  Inlet 

Best  Ramjet  * |  Best  IR  * 

Airframe  drag 
coefficient,  C 

Additive  drag 
coefficient,  C 


D 


0.349    0.349 


dad 


Thrust  coefficient,  Cf    0.349 


Combustor  exit  temperature 


R,  T 


t3 


4122 


Mass  flow  rate,  lb  /sec    9.5 

m 


67.2 


Specific  Thrust, 
lb£/(lbm/sec) 

Specific  fuel  consumption, 

(lb  /hr)/lb.  2.02 

m       t 

Maximum  thrust  coefficient 
C,-,  at  stoichiometric 
fuel  air  ratio  0.505 


Combustor  exit  tempera- 
ture, T  -,  °R 

Relative  detection 
range 


5698 


0.33 


0.018 
0.349 

3284 
18.3 

34.8 

2.68 

0.798 

5698 

0.29 


0.349 

0.085 
0.349 

4422 
13.6 

46.8 

3.35 

0.509 
5698 

0.37 


Best  ramjet  is  obtained  for  r,/r   equal  to  1.854  at  least 
for  the  range  of  rL/r   investigated  here.   Best  IR 
detection  capability  occurs  for  r,/rn  equal  1.63. 


n 
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V.   CONCLUSIONS 

Conversion  of  the  Navy  SALGP  from  a  rocket  to  a  ramjet 
severely  degrades  detection  capability,  if  an  axisymmetric 
inlet  is  used.   Note  the  values  of  0.33,  0.29,  and  0.37  for 
relative  detection  range  from  Table  IV-1. 

The  ramjets  with  blunt  nose  inlets  suffer  in  both  speci- 
fic thrust  and  specific  fuel  consumption.   The  poor  perfor- 
mance is  due  to  two  causes.   First,  the  high  value  of  C,  _, 
and  second,  poor  pressure  recovery  of  the  inlet  leading  to 
poor  performance. 

Surprisingly,  the  excess  thrust  coefficient  of  ramjets 
with  blunt  noses  is  very  competitive  with  the  ramjet  with 
spike  inlet. 

In  view  of  the  foregoing  conclusions,  pop-out  or  retrac- 
table scoop  inlets  appear  much  more  attractive  in  spite  of 
structural  constraints  and  space  considerations.   Use  of  an 
axisymmetric  inlet  located  in  the  nose  of  a  5-inch  projectile 
is  a  severe  design  compromise. 
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APPENDIX  A 
NASA  AMES  COMPUTER  PROGRAM  DESCRIPTION 

I.  INTRODUCTION 

Two  different  computer  programs  were  used  to  model  the 
flow  around  the  blunt  nosed  body.   The  coordinate  system 
used  in  the  thesis  and  the  NASA  Ames  computer  program  is 
shown  in  Figure  II- 6. 

The  first,  IMPLCBO,  is  a  modified  version  of  AXI -BLUNT; 
see  Kutler,  Chakravarthy ,  and  Lombard  [Ref.  9].   Programmed 
to  predict  the  supersonic  flow  over  a  three  dimensional 
body,  the  algorithm  predicts  shock  shape  and  location  as  well 
as  flow  parameters  at  equally  spaced  interior  points  between 
the  body  and  the  shock. 

The  second  program,  OGIVE,  solves  for  the  supersonic  flow 
parameters  around  a  three  dimensional  wing  body  configuration; 
refer  to  Kutler,  Reinhardt,  and  Warming  [Ref.  10].   Originally 
designed  for  flow  prediction  around  a  delta  winged  spacecraft, 
the  program  will  provide  flow  field  parameters  at  various  radii 
from  the  centerbody  as  well  as  at  various  locations  along 
the  body. 

II.  PROGRAM  DESCRIPTION 

The  IMPLCBO  program  uses  an  unsteady,  implicit  numerical 
procedure  to  determine  the  supersonic  flow  around  the  body. 
Further  program  description  is  contained  in  Reference  9. 
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The  program  output  is  referenced  to  various  free  stream 
values.   These  values  are  shown  in  Table  A-l. 

The  program  utilizes  the  body  configuration  and  predicted 
shock  points  as  boundary  values  and  solves  the  inviscid 
Euler  or  Navier-Stokes  equations  for  the  interior  points, 
after  a  coordinate  transform. 

The  different  angular  relationships  at  the  shock  front 
are  shown  in  Figure  A-l. 

The  OGIVE  program  utilizes  the  flow  field  parameters 
determined  by  IMPLCBO ,  at  9  =  90°  and  X  =  1.00,  as  an  initial 
starting  point  and  determines  flow  field  values  along  the 
cylinder  body  by  a  finite  difference  algorithm.   The  program 
is  further  described  in  Reference  10.   The  program  reference 
values  are  tabulated  in  Table  A-l. 

The  program  describes  the  body  shapes  by  analytical  approx- 
imations.   These  approximations  consist  of  coordinate  posi- 
tions and  the  slope  of  the  body  contour  as  a  function  of  the 
distance  along  the  body  axis. 
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Table  A-l 
Reference  Values  For  OGIVE  And  IMPLCBO  Programs 


IMPLCBO  (region  1) 

output  variable 

P/PINF 
RHO/RINF 
U/QINF 
V/QINF 


reference  value 


Poo 

^00 


ZyM    ;    at  M  =3.0 

— i-         00  '  00 

/7M    ;    at  M  =3.0 


OGIVE  (region  2) 

R 
P 


RHO 

QINF 
U 

V 

W 

Z 


r     =   1.0 
P?co 


q    =v   /v 

Hoo        oo'       m 


^00 

r 
n 


>    Poo/Ptoo=0-0272237    at 
Moo=3.0 

;  p  /p.  =0?0762263  at 


t°° 


M   =3.0 


•q0)  =  [l  +  5/Mi]1/2 

Z  component  of  flow 

velocity 

R  component  of  flow 

velocity 

0  component  of  flow 

velocity 

length  along 

cylinder  body 
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Bow  Shock 


Figure  A-l  Angular  Relationship  At  The 
Shock  Front  For  5,  S  ,  6 
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APPENDIX  B 
MASS  CONTINUITY  METHOD  OF  DETERMINING  STREAMLINE  SHAPE 

A  computer  program  has  been  written  for  the  HP  9830  to 
calculate  the  capture  streamline.   The  program  is  designed 
to  reduce  the  output  data  obtained  from  the  NASA  Ames  com- 
puter output  and  compute  the  streamtube  position  at  various 
points  in  the  flow  field  around  the  blunt  nosed  body. 

The  streamtube  position  is  determined  by  comparing  the 
mass  flow  of  air  entering  the  bow  shock  to  the  integral  of 

mass  flux.   The  streamtube  radius  at  the  bow  shock,  r  ,  is 

'   s 

determined  by  comparing  the  mass  flow  at  the  inlet  to  mass 
flow  at  the  nose. 

The  value  of  the  integral  in  equation  (22)  or  (23)  is 
determined  by  assuming  a  linear  change  in  pressure,  density 
and  velocity  between  successive  data  points.   The  program 
interpolates  to  find  the  distance  from  the  axis  that  the 
streamtube  must  be  to  satisfy  the  mass  flow  requirements  of 
the  inlet.   Values  of  X  and  Y  as  shown  in  Figure  1 1 -  6  are 
calculated,  and  an  interpolated  pressure  is  determined. 
The  value  of  the  integral  in  equation  (22)  and  (23)  is  cal- 
culated using  the  rectangular  method.   Equations  (22)  and 
(23)  are  developed  in  Chapter  III  of  this  thesis. 


74 


Table  B-l 
Variables  For  Capture  Streamtube  Computer  Program 

Symbol  for   Definition 
HP  9830 

AO  value  of  integral  in  equation  (22)  or  (23) 

between  point  i  and  i+1 

A7  remaining  area  between  tabulated  values 

C  squared  radius  of  streamtube  at  shock  front 

(square  inches) 

CI  constant  for  integral  when  in  region  I 

C2  constant  for  integral  when  in  region  II 

F  percentage  of  distance  Y7  is  between  Y(I) 

and  Y(I+1) 

G  dummy  variable  for  area  calculation 

I  indexing  variable 

P7  interpolated  value  of  pressure  X7,Y7 

S7  slope  of  area  curve 

W  dummy  variable  used  in  calculation  of  C 

Q9  value  of  Q^,  equal  to  (1+  (5/M^))"1^ 

R9  value  of  p  /p^   at  M 

S  summation  of  integral  area  to  the  ith  point 

(inches  squared) 

X7  interpolated  value  of  X,  determined  from  Y7 

(inches) 

Y7  interpolated  value  of  Y  (inches) 

Input  variables 

A  angular  value  of  data  points  above  projectile 

centerline  (degrees) 

19  number  of  data  points  in  data  file 

Z  location  of  data  points  in  relation  to  the 

projectile  shoulder.   1=  in  front  of  shoulder 
2=  on  the  body  cylinder 
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Table  B-l  con't 

Symbol  for.   Definition 
HP  983Q 

Dimensioned  Variables 

G(I)  integral  value  between  i  and  i+1 

L(I)  integral  value  at  ith  point 

P(I)  pressure  ratio  at  ith  point 

R(I)  density  ratio  at  ith  point 

U(I)         component  of  flow  velocity  along  the  body 
in  the  Z-direction 

V(I)         component  of  flow  velocity  along  the  body 
in  the  r  direction 

Y(I)         distance  from  the  body  centerline  to  the 
data  point 
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APPENDIX  C 
ADDITIVE  DRAG  COMPUTATION 

The  program  is  designed  to  calculate  the  additive  drag 
and  coefficient  of  additive  drag  for  an  inlet.   The  reference 
area  is  the  projectile  base  area.   The  program  is  written  in 
BASIC  for  the  HP  9830  desk  top  calculator. 

The  input  data  consists  of  the  point  position  (X  and  Y) 
and  the  static  pressure  at  the  point.   The  three  variables 
X,  Y,  and  p  are  determined  by  first  running  the  program 
described  in  Appendix  B  of  this  thesis. 

Equations  (27)  through  (30)  used  in  this  calculation  are 
developed  in  Section  III  of  this  thesis. 

The  computer  symbols  are  listed  and  defined  in  Table  C-l 
and  the  program  listing  is  in  Table  C-2. 
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Table  C-l 
Variables  For  Coefficient  Of  Additive  Drag  Computer  Program 

Variables    Definition 

A  Reference  area  (square  inches) 

C  Coefficient  of  additive  drag  between  I  and  1+1 

data  point  (dimensionless) 

CI  total  coefficient  of  additive  drag  (dimensionless) 

D  total  drag  of  body 

H  distance  along  body,  Z,  between  data  points 

(inches) 

I  counting  variable 

J  counting  variable 

P  Ratio  of  circle  circumference  to  diameter 

PI  p   /p„      at   input  M 

^OO'     r  £00  C  00 

R  Distance  from  centerline,  Y  (inches) 

S  Slant  distance  between  data  points  (inches) 

SI  Surface  area  of  cone  between  data  points 

(square  inches) 

W  Incremental  drag  between  data  points  (lb) 

Input  variables 

M  Mach  number  (dimensionless) 

19  Number  of  data  points 

Dimensioned  variables 

A(I)         Angle  of  streamtube  between  data  points  (degrees) 

P(I)         Pressure  at  Ith  data  point  (dimensionless) 

R(I)         Average  pressure  between  data  points 
(dimensionless) 

X(I)         X  position  of  Ith  data  point  (inches) 

Y(I)         Y  position  of  Ith  data  point  (inches) 
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